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ABSTRACT 



Context. The spectacular outflow of HH 46/47 is driven by HH 46 IRS 1, an embedded Class I Young Stellar Object (YSO). Although 
much is known about this region from extensive optical and infrared observations, the properties of its protostellar envelope and 
molecular outflow are poorly constrained. 

Aims. Our aim is to characterize the size, mass, density and temperature profiles of the protostellar envelope of HH 46 IRS 1 and its 
surrounding cloud material as well as the effect the outflow has on its environment. 

Methods. The newly commisioned CHAMP + and LABOCA arrays on the APEX telescope, combined with lower frequency line 
receivers, are used to obtain a large (5' x 5', 0.6x0.6 pc) continuum map and smaller (80" x 80", 36,000x36,000 AU) heterodyne 
maps in various isotopologues of CO and HCO + . The high-./ lines of CO (6-5 and 7-6) and its isotopologues together with [C I] 2-1, 
observed with CHAMP + , are used to probe the warm molecular gas in the inner few hundred AU and in the outflowing gas. The data 
are interpreted with continuum and line radiative transfer models. 

Results. Broad outflow wings are seen in CO low- and high-/ lines at several positions, constraining the gas temperatures to a 
constant value of ~100 K along the red outflow axis and to ~60 K for the blue outflow. The derived outflow mass is of order 0.4—0.8 
M Q , significantly higher than previously found. The bulk of the strong high-/ CO line emission has a surprisingly narrow width, 
however, even at outflow positions. These lines cannot be fit by a passively heated model of the HH 46 IRS envelope. We propose that 
it originates from photon heating of the outflow cavity walls by ultraviolet photons originating in outflow shocks and the accretion 
disk boundary layers. At the position of the bow shock itself, the UV photons are energetic enough to dissociate CO. The envelope 
mass of ~5 M Q is strongly concentrated towards HH 46 IRS with a density power law of -1.8. 

Conclusions. The fast mapping speed offered by CHAMP + allows the use of high-/ CO lines and their isotopes to generate new 
insights into the physics of the interplay between the molecular outflow and protostellar envelope around low-mass protostars. The 
UV radiation inferred from the high-J CO and [C I] data will affect the chemistry of other species. 



1. Introduction 

The Young Stellar Object (YSO) HH 46 IRS 1 (RA= 
08h25m43.9s, Dec =-51d00h36s (J2000) ), located at the edge 
of an isolat ed Bok Globule (D = 450 pc) in the southern 
hemisphere dSchwartzl 1 1977]). is well-k nown for its spectacu- 
lar outflow (jReipurth & Heat hcotdl 19911) . observed at both visi- 
ble and infrared (IR) wav elengths with the Hubb l e and Spitzer 
Space Telescopes (e.g.. iHeathcote et alj 119961 IStanke et al.l 



1999, iNoriega-Crespo et alj|2004 IVelusamv etall 120071) . Deep 
Ho- observations using the VLT have revealed bow shocks as- 
sociated with t he HH 46 outflow up to a parsec away from the 
central source dStanke et al.lll999 l). Its blue-shifted lobe expands 
outside the cloud in a low density region, due to the close prox- 
imity of the protostar to the edge of the cloud. Proper motion 
and radial velocity studies show that the outflow has an inclina- 
tion of 35° with respect to the plane of th e sky and flow ve- 
locities in atomic lines up to 300 km s _1 dDopita et al.l fl982. 
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Reipurth & HeathcoBQlH iMicono et alJ Tl998>. The internal 
driving source was fo und to be HH 46 IRS 1 (L = 16 L ), an em- 
bedde d Class I YSO dRaymond et alJll994ISchwartz & Greene! 
2003). Surprisingly, much less is known about the proper- 
ties of the protostellar e nvel ope and the mole cular outflow. 
IChernin & Massonl d!991l) and lOlberget all d 19921) mapped this 
region using low excitation CO lines, which show that contrary 
to the optical flows, the red -shifted outflow lobe is muc h stronger 
than the blue-shifted one. IChernin & Masson] d!991l) theorised 
that a lack of dense material in the path of the blue-shifted flow 
is responsible for this. 

Comparisons between dust emission and molecular lines at 
submillimeter wavelengths, together with self-consistent radia- 
tive transfer calculations, have been extensively used to charac- 
terize the phy sical and chemical structure of Class and Class 
I envelopes dSchoieret all 120021 LMreensen etall 120021; l2005bl 
iMaret et alJ,l2004l) . However, an essential component could not 
be probed with those data. The amount of warm (T > 50 K) gas 
within the protostellar envelope as well as the influence of the 
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molecular outflow have not been constrained directly using ob- 
servations of lower excited molecular lines in the 200 and 300 
GHz bands. Although more complex molecules, such as H2CO 
and CH3OH emit at these frequencies from high en ergy levels 
(e.g. [van Dish oeck et al.ll995llCeccarelh et al.l20 00). their more 
complex chemistry complicates their use as probes of the warm 
gas. Observations of CO at higher energies (up to 200 K) pro- 
vide more reliable probes into the inner regions of envelopes and 
molecular outflows, but su ch lines have only been observed for 
a han dful of sources (e.g., Hjo gerheiid eet al.lll998l iParise et alj 
[20061 1 van Kempen et al.l l 2006). [C I] emission provides an im- 
portant co nstraint on the stren gth of the radiation field within the 
outflows dWalker et al.lll993b . 

With the commi ssioning of the Atacama Pathfinder 
Expe riment (APEX j^ Giiste n et al.ll2006l). the CH AMP+ instru- 
ment (Kasem ann et alj|2006l iGiisten et alj|2008l) allows obser- 
vations of molecular emission lines in the higher frequency sub- 
millimeter bands of southern sources, like HH 46. CHAMP + is 
the first array of its kind. With its 14 pixels, it is able to ob- 
serve simultaneously in the 690 and 800 GHz atmospheric win- 
dows. This combination of dual-frequency observing and fast 
mapping speed, supplemented by lower fre quency single pixel 
data and a LA BOCA continuum array map dKrevsa et al.| [2003. 
Sirin go et al.l 120081) . provides a large range of highly comple- 
mentary tracers of both the gas and dust conditions in the in- 
ner and outer regions of the envelope, as well as the molecu- 
lar outflow on scales of a few arcminutes. Spectral line maps 
provide key information that is essential in the analysis of em- 
bedded YSOs which single-pointed observa tions cannot offer 
dBoogert et alj|2002llvan Kempen et af]l2008l) . 

In this paper, we present first results from CHAMP + and 
LABOCA observations of the HH 46 IRS 1 source, supple- 
mented by lower frequency line receivers. Observations of 
highly excited CO, HCO + and their isotopologues are used to 
constrain the properties of the protostellar envelope and molecu- 
lar outflow of HH 46. In § 2 the observations performed at APEX 
are presented. Resulting spectra and maps can be found in § 3. 
In § 4 and 5 the protostellar envelope, molecular outflow and 
close surrounding of HH 46 are characterized using a radiative 
transfer analysis. In addition, we discuss a possible scenario for 
the high-/ CO emission in § 6. The final conclusions are given 
in § 7. CHAMP + observations of a larger sample of low-mass 
protostars are presented in a subsequent paper (van Kempen et 
al. submitted). 



2. Observations 

Mole cular line observations were carried out with the CHAMP + 
array dKasemann et al. 2006) of CO and its isotopologues, rang- 
ing in transitions from 7=6-5 to 7=8-7, as well as [C I] 2-1. 
These observations were supplemented with low excitation line 
observations using APEX-1 (230 GHz, CO 2-1), APEX-2a (345 
GHz, CO, C ls O and 13 CO 3-2, and HCO + and H 13 CO + 4-3) and 
FLASH (460 GHz, CO and 13 CO 4-3). In addition, LABOCA 
was used to map the entire region at 870 /mi. See Table Q] for 
an overview of the observed emission lines for each instrument 
and the corresponding rest frequencies and upper level ener- 



1 This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX) in programs E-77.C-0217, X-77.C- 
0003, X-79.C-0101 and E-081.F-9837A. APEX is a collaboration 
between the Max-Planck-Institut fur Radioastronomie, the European 
Southern Observatory, and the Onsala Space Observatory. APEX-1 was 
used during science verification in June 2008. 



gies, together with their corresponding APEX beams. All ob- 
servations were done under excellent weather conditions with 
typical system temperatures of 2100 K for CHAMP+-I (SSB, 
690 GHz), 7500 K for CHAMP+-II (SSB, 800 GHz), 1100 K 
for FLASH-I (DSB, 460 GHz), 290 K for APEX-1 and 230 
K for APEX-2a (both SSB). Calibration errors are estimatd at 
15 to 20 %. The HH 46 protostar was spectrally mapped in 
CO 2-1, 3-2, 4-3, 6-5, 7-6 and HCO + 4-3, as weU as 13 CO 
6-5 and [C I] 2-1. The mapped area differs per line, ranging 
from 40" x 40" for CO 4-3 to 200" x 200" for CO 2-1, with 
most other lines covering 80" x 80". Observations were taken 
over a period of 2 years from July 2006 to Septe mber 2008 
using Fast Fourier Transform Spectrometer (FFTS) dKlein et al.l 
2006) backends for all instruments, except CHAMP + , for which 
only the two central pixels were attached to the FFTS back- 
ends. Other CHAMP + pixels were attached to the MPI Array 
Correlator System (MACS) backends. FFTS backends are able 
to reach resolutions of 0.12 MHz (0.045 km s 1 at 800 GHz), 
while the MACS units were used at a resolution of 1 MHz (0.36 
km s _1 at 800 GHz). Beam efficiencies are 0.75 for APEX-1, 
0.73 for APEX-2a, 0.7 for FLASH-I, 0.56 for CHAMP+-I and 
0.43 for CHAMP + -II. Pointing was checked on nearby sources 
and found to be accurate within 3" for the APEX-2a observa- 
tions. For CHAMP + , pointing is accurate within ~5". All ob- 
servations were taken using position switching with reference 
positions in azimuth ranging from 600" to 3600". 

LABOCA observed HH 46 during November 2007 using its 
295 pixels in a spiral mode centered on the IR position using a 
position switch of 600" arcseconds. Only the inner 5' x 5' of 
the 11' field of view was used. The field was integrated down to 
a noise level of 0.01 Jy/beam, averaged over the entire region. 
Due to the used spiral mode, the map contains a radial increase 
of noise towards the edge of the map. The continuum data were 
reduced with the BOA packag^l 

3. Results 

3.1. Single pixel spectra 

Figure Q] shows the spectra taken at the position of HH 46 IRS 
1 for all lines in Table Q] Integrated intensities, peak tempera- 
tures and estimated contributions from the blue- and red-shifted 
outflowing gas are given in Table [2] The latter are derived by 
subtracting the central part of the line profile, associated with 
the quiescent gas, with a single gaussian. Emission was detected 
for all lines with the exception of C ls O 6-5 and 13 CO 8-7. The 
quiescent gas component peaks at a Vlsr of 5.3 +0.1 km s and 
has a FWHM of 1-2 km s _1 depending on the energy of the up- 
per level. Since the detection of the [C I] 2-1 line is only 3.5 <x, 
a gaussian fit is overplotted in Fig.Q] 

Integrated intensities range from 82.5 K km s _1 for CO 3-2 
to 1.1 K km s _I for H 13 CO + . All 12 CO line profiles show contri- 
butions of a red-shifted outflow lobe within the beam, including 
the high excitation CO 7-6 line. Emission from the blue-shifted 
outflow lobe is much weaker and not found for CO 7-6. In the 
other lines, outflow emission is only detected for HCO + 4-3, 
where a weak red-shifted wing is found. 

3.2. Maps around IRS 1 

The resulting dust emission map (5' x 5') at 870 //m is shown in 
Fig. [2] The envelope is clearly resolved, considering the beam- 

2 http://www.astro.uni-bonn.de/boawiki 
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CO 
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16.6 
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28 
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200" X 170" 




3-2 


33.2 


345.796 


APEX-2a 
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80" X 80" 




4-3 


55.3 


461.041 


FLASH-I 


14 
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40" X 40" 




6-5 


116.2 


691.473 


CHAMP+-I 
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80" x 80"° 




7-6 


154.8 


806.651 


CHAMP + -II 
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80" x 80" 
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n 
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50" x 50" 
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CHAMP" 1 " -II 
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n 




c 18 o 
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31.6 
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APEX-2a 


18 


n 






6-5 


110.6 


658.553 


CHAMP+-I 


9 


n 




[C I] 


2-1 


62.3 


809.342 


CHAMP+-II 


8 


y 


50" x 50"" 


HCO+ 


4-3 


42.8 


356.724 


APEX-2a 


18 


y 


80" x 80" 


H 13 CO+ 


4-3 


41.6 


346.998 


APEX-2a 


18 


n 






870 pm continuum 




LABOCA 


18 


y 


6' x6' 



" additional mapping was performed in these lines around the bow shock position. 
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Fig. 1. Single spectra taken at the central position of HH 46 (all in order from bottom to top). Left : 12 CO 2-1, 12 CO 3-2, 12 CO 4-3, 
12 CO 6-5 and 12 CO 7-6. Middle : 13 CO 3-2, 13 CO 4-3, I3 CO 6-5 and I3 CO 8-7. Right : C I8 3-2, C I8 6-5, H I3 CO + 4-3, HCO + 
4-3 and [C I] 2-1. Spectra have been shifted vertically for easy viewing. 



HH 46 CO 3-3 HH 46 c0 7 _ 6 HH46 HCO* 4-3 




A RA ["] A RA ["] A RA ["] 

Fig. 3. From left to right : Spectrally integrated CO 3-2, 7-6 and HCO + 4-3 maps of HH 46 (see Table|2|. Contour levels are at 3cr, 
6cr, 9cr, ... (cr= 1.5 K for CO 3-2, 0.4 K for HCO + 4-3 and 0.8 K for CO 7-6) The position of the IR source is marked with a X. 
Note that the HCO + 4-3 map has higher noise levels than the central spectrum in Fig[TJ 

size of 18" of the APEX dish. The envelope is slightly elongated on a south-west to north-east axis. South of HH 46, dust emis- 
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Table 2. Observed molecular line intenstities". 
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CO 


2-1 


53.3 


22.4 


11.5 


3.8 


61.7 


21.2 




3-2 


82.5 


19.4 


28.0 


13.4 


65.6 


19.8 




4—3 


70.5 


14.9 


44.2 


3.7 


38.5 


14.0 




6-5 


42.6 


9.5 


11.7 


5.7 


34.6 


8.8 




n t 
l—O 


40. J 


8.6 


17.6 


<i. I 


25.3 


5.9 


CO 


3-2 


14.3 


11.4 












4-3 


8.0 


5.5 












6-5 


3.9 


2.9 






<0.6 rf 
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3.2 


3.3 












6-5 




<0.5 rf 










[CI] 


2-1 


2.3 


2.7 








<l.S d 


HCO+ 


4-3 


8.5 


5.4 


1.4 








H 13 CO + 


4-3 


1.1 


1.0 











" Red and Blue shifted emission are calculated by subtracting the central part of the line profile with a gaussian with FWHM of 1.5-2 km s 

6 Position : ARA = -20", ADec = -20". 

Velocity range used for integration: -5 km s _1 to 15 km s -1 

d Upper limit of 3<x in a 0.7 km s _I channel. 



M Q , while the central 120" is associated with a mass of 3.7 M , 
both with an assumed dust temperature of 20 K. 

The integrated spectral line maps of CO 3-2, 7-6 and HCO + 
4-3 (all 80" x 80") in Fig.[3]show a similar structure as the con- 
tinuum maps (see also §5.4.1), although the elongation to the 
south-west is much more pronounced in the CO lines, with an 
integrated intensity of CO 3-2 of 65.6 K km s _1 seen in the - 
20",-20" position. As can be seen from the red- and blue-shifted 
velocity maps of CO 2-1, 4-3 and 6-5 in Fig. H] the shape of 
the integrated intensity maps is largely due to the outflow con- 
tributions, especially for the low-excitation lines, even in the CO 
2-1 map, for which a map of 200 " x 200" was observed. Th is is 
in agreement with the results from lChernin & Masson (119911) . In 
fact, the CO 3-2 line is so prominent that it can contribute sig- 
nificantly to the 'continuum' emission seen by LABOCA. The 
CO 3-2 line produces a flux density of 40 mJy/beam off source, 
equivalent to the LABOCA 4<x level. The total continuum emis- 
sion seen at these positions is ~12-15<x. CO emission may thus 
contribute up to 30% to the observed LABOCA emission in the 
outflow region. Another possible contribution to the dust emis- 
sion at 870 /mi is heating of the dust grains by UV radiation that 
takes place in the cavity walls. This is further discussed in § 6. 

In contrast, the elongation in the HCO + 4-3 map is much 
less pronounced than in the CO low-/ data. The map of [C I] 
2-1 is not shown, since no lines were detected down to 1.8 K 
(3(f), except at the source position (Fig. 1) and at the bow shock 
(see § 3.3). In the map of 13 CO 6-5, lines are only detected at the 
central position and at neighbouring pixels along the outflow. 

Fig.|5]shows the total distribution of 12 CO 6-5 over the entire 
area, o verplotted over the Spitzer image of Noriega-Crespo ^FaT] 
(120041) with the Spitzer-IRAC 1,2 and 4 bands (3.6, 4.5 and 8.0 
fim respectively.) The CO 6-5 integrated intensities follow the 
outflow but do not extend all the way out to the bow shock. In 
addition, several maximum intensities can be found. The first 
one corresponds with the protostellar envelope, but the other 
two seem to be related to the outflow. Similarly regularly spaced 
'knots' are seen in , for example, the CO 3-2 maps of the NGC 
1333 IRAS4A outflow (BlakeetaL|[i^l. 




150 100 



50 -50 

RA offset ["] 



-100 -150 



Fig. 2. Continuum map at 870 /im of HH 46 obtained with 
LABOCA with contour levels plotted in multiples of the 3<x 
noise level of 0.05 Jy/beam. Left inset: the radial profile going 
outwards from the IR position. Right inset: the SED from HH 46 
from 10 /im to 10 mm. In both panels, the best-fitting DUSTY 
model is overplotted. 



sion from the cold cloud is seen. The total integrated flux in a 
120" diameter aperture around the source position is found to 
be 3.7 Jy. Reduction of archiv al SCUBA (850 um) data of this 
region yielded a flux of 3.3 Jy dDi Francesco et ai1l2008l) . a dif- 
ference in flux well within the estimated calibrat ion errors of 
both S CUBA and LABOCA. Using the formula in lShirlev etafl 
(2000), the total mass in the 300" mapped area of the cloud is 8 
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HH 46 CO 2-1 HH 46 CO 6-5 HH 46 CO 4-3 




4 RA ["] A RA ["] A RA ["] 

Fig. 4. Maps of HH 46 showing the outflowing gas. Solid lines show blue-shifted, dashed lines how red-shifted emission. Contours 
are drawn at 3cr, 6<x, 9cr,... , with the exception of CO 2-1, where contours are drawn 3<x, 9<x, 15cr, ... . 3cr levels are 0.4 K km s , 
1 K km s" 1 and 1 K km s _1 for 2-1, 6-5 and 4-3 respectively. Blue- and red-shifted emission is derived by subtracting a central 
gaussian fitted to the quiescent central part of the line profiles. Velocities typically range from -5 to 1.5 km s for the blue-shifted 
emission and 8 to 15 km s _1 for the red-shifted. Note that the images are at different spatial scales and progressively zoom in onto 
the central protostar from left to right. The position of the infrared source is marked in all maps with a 'X'. 




Fig. 5. The integrated intensity of CO 6-5 (blue contours) and [C I] 2-1 (green contours) overplotted on the Spitzer-IRAC 1 (blue), 
2 (green) and 4 (red) bands of the entire HH46 region. CO 6-5 contours are in increasing order of 5 K km s _1 . The clear cut between 
the CO and [C I] emission near the bow shock suggests that UV photons capable of dissociating CO are present near the bow shock. 



Figure|6]shows the observed spectra within the central 80" x 
80" of CO 3-2, taken with APEX-2a, and CO 6-5, taken with 
CHAMP+-I. Both are binned to square 10" x 10" pixels. The line 
profiles show an interesting distribution of emission, especially 
at the positions associated with the outflowing gas. The CO 3- 
2 spectra consist of a central gaussian originating in the cold 
envelope material surrounding the protostar, with a strong red- 
shifted outflow of up to 10 km s away from the source and 



cloud velocity in the south-east direction. The line profile at the 
source position has a quiescent component, flanked by both red- 
and blue-shifted emission. Interestingly, the CO 6-5 emission 
shows a quiescent component with relatively weak red-shifted 
emission at other positions in the map. The central part of the 
line profile can be fitted with narrow (AV = 2-3 km s~') gaussian 
profiles. At the central position, outflow emission in the CO 6-5 
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HH 46 CO 3-2 HH 46 CO 6-5 




Fig. 6. Spectra of HH 46 of CO 3-2 (left) and CO 6-5 {right) in the 80" x 80" mapping area. Individual spectra are shown on a 7mb 
scale of -5 to 20 K for CO 3-2 and from -3 to 12 K for CO 6-5. Both axes are from -10 to 20 km s _1 . The black lines show the axis 
of the red outflow. 



line is much more prominent, but a strong quiescent component 
is still present. 




-10 10 20 

Velocity (km/s) 



Fig. 7. Single spectra of (from bottom to top): CO 2-1, CO 3-2, 
HCO + 4-3, CO 4-3, I3 CO 6-5, CO 6-5, [C I] 2-1 and CO 7-6. 
All spectra are at an position in the red outflow lobe of ARA=- 
20" and ADec=-20". The quiescent gas has a FWHM of 1.5 km 
s _1 centered at 5.3 km s _1 . 

Figure [7] shows all spectra observed within the spectral line 
maps at a relative position of (-20", -20"). Peak temperatures 
and integrated intensities are given in Table|2] This position cov- 
ers the red-shifted outflow seen prominently in the emission of 



the low-excitation CO 3-2 and 2-1 lines. Even CO 4-3 shows 
significant red-shifted outflow emission. The outflowing gas is 
still present in the high-/ CO lines. However, for both CO 6-5 
and 7-6 a significant part of the emission (on the order of 70- 
80%) originates in a quiescent narrow component. The isotopo- 
logue 13 CO 6-5 is not detected, down to a 3<x level of 0.6 K in a 
0.7 km s _1 bin. Similarly, no [C I] is detected down to 1.8 K and 
no HCO + 4-3 is detected down to 0.9 K, both 3<x in a 0.7 km s _1 
bin. 

Using the limits on the 13 CO 6-5 emission at both the cen- 
tral position and the selected off-postion, it is found that the 
quiescent component at the central position is optically thick 
(t > 3), while the quiescent emission at the outflow position 
is optically thin (t < 0.4), as no 13 CO is detected between 
and 10 km s" 1 . This an alysis assumes a 12 CO: 13 CO ratio of 70:1 
dWilson&Roodlll994l) . 

Although outflow emission heavily influences the line pro- 
files of the CO 3-2 and 6-5 throughout the maps, positions 
south, south-east, and north of the (0,0) position are not affected 
by any outflowing gas, as seen in Fig. [6] About 30" to 40" north 
of HH 46, both the CO 3-2 and CO 6-5 are not detected. Even 
in the map of CO 2-1, no emission was found at these positions. 
It is concluded that the cold cloud material of the surrounding 
Bok globule does not extend to these scales. At positions south 
and south-east of HH 46, CO 3-2 emission is seen, but no CO 
6-5 is detected there. 

3.3. Bow shock 

The bow shock associate d with the red outflow lobe i s clearly 
visible in the IR images of Norieg a^Crespo et"ai1 (120041) . located 
at a relative offset of (-100", -60") with respect to the source. 
Additional observations of this bow shock in CO 6-5, [C I] 2- 
1 and 13 CO 6-5 were carried out over a 90" by 60" area. Fig. 
[5] shows the distribution of [C I] 2-1 emission in this region. 
Narrow [C I] 2-1 is clearly detected near and at the position of 
the bow shock, where it has a integrated intensity of 1.2 K km 
s _1 with a peak temperature of 1.7 K. Interestingly, the emission 
seems to be spatially extended along the outflow axis into the 
outflow lobe. No I3 CO 6-5 was detected at any position near the 
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bow shock down to an rms of 85 mK in a 0.3 km s _1 bin. The 
combination of detected [C I] 2-1 and lack of CO 6-5 emission 
suggests that CO is dissociated by either the shocks present near 
the bow shock or by UV photons capable of dissociating CO )see 
§6.1). 



4. Envelope and surrounding cloud 

4. 1 . Envelope - dust 
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Fig. 8. The density (left) and temperature (right) distributions of 
the DUSTY modelling of HH 46. 



Using the 1-D radiative transfer code DUSTY 
dlvezic & Elitzurl [l997), a spherically symmetric envelope 
model is constructed by fitting the radial profile of the 870 fim 
image and the SED simultaneously, determining the size, total 
mass, inner radius as well as density and temperature profiles 
of the protostellar e nvelope. For a more t horough discussion 
of this method, see |j0rgensen etaH <[2002). DUSTY uses / 
(=^outer/^innerX P> me power law exponent of the density 
gradient (n oc r~ p ), and tioo, the opacity at 100 fim as its free 
parameters. The temperature at the inner boundary was taken 
to be 250 K. For the SEP flux es were used at 60 and 100 
urn (IRAS. [Henning et al ] [1993b, 850 fim (SCUBA archive, 
iDi Francesco et al.l I2008L 870 um ( LABOCA, this work) and 



1.3 mm (SEST lHenning etal.1119931). No emission from MIPS 
at 24 fim or ISO-SWS dNisini et alJl20'ol iNoriega-Crespo et alJ 



2004) was used. The 24 fim flux is included in the figure as 
reference. Note that the model is unable to account for this high 
flux. Deviations from sp herical symmetry (ou tflow cavities), 
additi onal shock emission (fVe lusamv et al. 2007) or larger inner 
holes dj0rgensen et alj|2005ah are often inferred to explain the 
observed high mid-IR fluxes that cannot be fitted with DUSTY. 
The radial profile was determined in directions away from the 
southwest outflow and cloud material, ignoring any emission 
in a 90° cone to the south-west. Results for the best-fitting 
envelope model to the 16 L for the source luminosity can be 
found in Table [3] Temperature and density distributions are 
displayed in Figure [8] The corresponding fits are shown in the 
insets of Fig . |2] 

The envelope contains a large amount of cold gas (~5.1 M ) 
within the outer radius of 20,800 AU (~0.1 pc), but with a sig- 
nificant fraction of the mass concentrated towards the inner en- 
velope due to the steep density profile (p=1.8). The H2 density 
at 1,000 AU is 2.4xl0 5 cm 4 . The transition from envelope to 
parental cloud is not taken into account. 



Table 3. Results from radiative transfer modelling for envelope 
properties 



Best fit 



Envelope Properties 



Y 


700 


iW(250 K) 


30 AU 


P 


1.8 


fi(30 K) 


533 AU 


T"100 


2.6 


R(n ds ) 


6,100 AU 








20,800 AU 






/i(1000 AU) 


2.5X10 6 cm" 3 






M env (> 10 K) 


5.1 M 



4.2. Envelope - gas 

The physical structure of the gas is best traced by optically thin 
emission that probes the quiescent envelope gas at high density 
dj0rgensen et al.ll2005bl) . To that end, the temperature and den- 
sity structure derived from the dust radiative transfer model was 
used as a model in put for the model using the RATR AN radia- 
tive transfer code dHogerheiide & van der Tak 2000) with data 
files from the LAMDA database dSchoier et al.l l2005). Two dif- 
ferent scenarios were investigated . In the first, only freeze-out is 
taken into account, with an inner CO abundance Xo of 2.7xl0 -4 
with respect to H2 and an outer CO abundance Xj of 10~ 5 . In 
the second , the abundances are the same, except that a high 
abundance Xo=2.7xlO~ 4 has been adopted for the outermost en- 
velope regions where density is lower than «d e =10 5 cirT 3 . This 
so-called 'drop' abundance profile is motivated by the fact that 
at low densitites, the timescales for freeze-out onto the grains 
are longer than the typical lifetimes of the cor es of a few 10 s 
yr. Th e ab undances Xq and Xq wer e derived by J0rg ensen et al.l 
(2002) and J0rgens en et al.l d2005bl) for a range of sources, based 
on e mission from optically th in lines. Isotope ratios were taken 
from IWilson & Roodl <fl99l of 550 for CO:C ls O and 70 for 
CO: 13 CO. The velocity field of the envelope is represented by 
a turbulent width of 0.8 km s . The gas and dust temperature 
were assumed to couple throughout the envelope. Gas and dust 
temperatures decouple in the cold outer envelope region with 
the gas temperature dropping by a factor two with respect to th e 
dust temperature dCeccarelli et al.lll996llDotv & Neufeldll 19971) . 
This drop can be counteracted, however, by UV radiation from 
the outside impinging on the envelope and heating the gas. For a 
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more thorough discussion, see also J0rgens e~n et al.l d2002l) . Any 
small temperature difference in the outer envelope affecst mostly 
the low-/ 12 CO lines. 

Table 2] gives the results for these models. For each model, 
the envelope contribution of both J Tmb and r pea k are given. The 
optically thin lines (C ls O 3-2 and 6-5, as well as the B CO 8-7) 
show that the inclusion of a higher abundance outer layer, as in 
a drop model, is necessary to increase the C 18 3-2 emission to 
the observed levels. At the same time this does not change the 
emission of the C 18 6-5 line. The adopted abundances of the 
drop abundance (Model 2) agree best with the observed intensi- 
ties of the optically thin lines, including the upper limits on C ls O 
6-5 and 13 CO 8-7 within the uncertainties. 



Table 4. Results from molecular line radiative transfer mod- 
elling. 



Line 


Model 1° 


Model 2° 


Obs. 




f Tmb 


y b 
l peak 


f Tmb 


y b 
1 peak 


^peak 


CO 2-1 


9.4 


5.0/22.2 


10.4 


4.5/78.9 


22.4 


CO 3-2 


8.0 


4.7/22.5 


7.1 


3.8/55.8 


19.4 


CO 4-3 


7.3 


4.2/16.9 


4.7 


3.1/47.3 


14.9 


CO 6-5 


6.8 


3.4/4.8 


5.1 


2.7/5.2 


9.5 


CO 7-6 


4.1 


2.2 


4.0 


2.2 


8.6 


"CO 3-2 


4.6 


3.7/3.8 


5.9 


3.6/10.7 


11.4 


13 CO 4-3 


3.5 


2.8 


4.2 


2.6/3.8 


5.5 


13 CO 6-5 


1.9 


1.2 


1.9 


1.2 


2.9 


13 CO 8-7 


0.8 


0.5 


0.8 


0.5 


<0.9 


C 18 3-2 


1.3 


1.2 


3.5 


2.9 


3.3 


C 18 6-5 


1.0 


0.7 


1.0 


0.7 


<0.5 


HCO +1 4-3 


3.5 


2.3/11.8 


2.8 


2.5/8.3 


5.4 


H 13 CO+ 4-3 


0.75 


0.7 


0.75 


0.7 


1.0 



"Model 1 has a jump abundance profile with X /X d of 2.7x10 4 /lxl0 5 . 
Model 2 has a drop abundance profile with X /X d /X of 
2.7xl0- 4 /ixl0- 5 /2.7xl0- 4 . 

* For some lines, the peak temperature is given as X/Y. The first value 
is the actual T pca ^ in the model profile. The second value refers to the 
peak of a gaussian fitted to the line wings of model profiles which show 
self-absorption. If only a single number is given, the modelled line is 
gaussian in nature and does not have self-absorption. 
c HCO + abundances of X and X d are 2.0xl0~ 8 and 3.0X10" 9 . 

The very optically thick low-/ I2 CO (2-1, 3-2 and 4-3) lines 
cannot be fitted by either model due to the lack of self-absorption 
in the observed lines, which is strongly present in the modelled 
profiles. To obtain a rough correction for self-absorption, gaus- 
sian fits were made to the line wings of CO, similar to the one 
in van Kempen et al. submitted.. Such gaussian fits provide only 
upper limits to the peak emission of CO, as the true CO emission 
is best fitted with an infall velocity (Schoier et al. 2002). 

Even with gaussian fits, the high-/ (6-5 and higher) lines 
of both 13 CO and 12 CO are severly underproduced by almost a 
factor 3. Note that their model line profiles do not show self- 
absorption. The emission in these lines is almost identical in the 
drop or jump abundance models. Increasing or decreasing the 
CO abundances Xq and X4 is not possible since the emission of 
optically thin lines such as C ls O 6-5 and 3-2 would then either 
be over- or underestimated. The main conclusion from the en- 
velope models is therefore that an additional, relatively optically 
thin but hot component unobscured by the warm envelope re- 
gion is needed to account for this emission. The origin of this 
hot component producing narrow highly excited CO line emis- 
sion will be discussed in § 5.6. 



The integrated [C I] 2-1 intensity (2.3 K km s~') can be re- 
produced with a constant abundance ratio C/H2 of 3-5xl0~ 7 (or 
about 0.1-3% of CO, depending on radius), typical for the dens- 
est molecular clouds. The [C I] line is optically thin, even for 
much higher abundances. A C abundance as high as 5xl(T 6 can 
be maintained by photodissociation of CO due to cos mic-ray in- 
duced UV photons deep inside the envelope (e.g. jFlower et al.l 
11994 . 

HCO + 4-3 shows little to no difference between Model 1 and 
2, because most of the emission traces gas denser than 10 5 ctrT 3 . 
There is a significant difference in /peak from the line profile vs. 
gaussian fits for main isotope line because of self-absorption. 
The optically thin H 13 CO + is slightly underproduced, similar to 
the C ls O 3-2 peak temperature. 

4.3. Surrounding cloud material 

The envelope of HH 46 is surrounded by cold, quiescent cloud 
material as evidenced by the LABOCA map and by the CO lines 
to the south-east (e.g., 20", -40"). If the cloud is assumed to 
be isothermal and homogeneous, the CO emission maps from 

2- 1 to 7-6 can be used to constrain t he its properties with th e 
radiative transfer program RADEX ( Ivan der Tak et al.l l2007r 3 l 
RADEX calculates the non-LTE excitation and line emission 
of molecules for a given temperature and density using an es- 
cape probability formulation for the radiative transfer. For gas at 
a constant temperature and density, such as in the surrounding 
cloud, the results from RADEX and RATRAN are comparable, 
but RADEX is easier and faster to use. A line width of 1.2 km 
s is used from the gaussian fits to the CO 3-2 emission at the 
cloud positions. The ratios of CO 3-2/6-5 peak temperatures at 
these cold cloud positions are at least 20, since CO 6-5 is not 
detected. RADEX simulations show that the CO 6-5 line has 
to be sub-thermally excited, which can be done when the gas is 
at low densities or at very low (/ <10 K) kinetic temperatures. 
The line ratios are best fitted by a cloud with a temperature of 
~14 K, a low density of a few times 10 3 cm -3 and a CO col- 
umn density of 10 17 -10 18 crrT 2 . The spatial distribution of the 
CO 3-2 and 6-5 lines also indicate that these cloud conditions 
extend homogeneously to at least 40" south of HH 46, the extent 
of both CO 3-2 and 6-5 map s. Maps of the op t ically thin C ls O 

3- 2, as have been done by Ivan Kempen et al.l (120081) for other 
sources, are needed to further constrain the column density and 
spatial structure of the surrounding cloud. 



5. Outflow 

5.1. Outflow temperature 

The spectral line maps clearly reveal the red-shifted outflow to 
the south-west and the blue-shifted outflow to the north-east. The 
large diff erence in extent betw een the two outflow lobes, already 
noted by Olbe rg etaD (Q992), is seen in all transitions, with the 
blue-shifted lobe producing much weaker emission. The brighter 
and larger outflow lobe seen with Spitzer dNoriega-Crespo et al.l 
2004, IVelusamy et al.ll2007l) corresponds to the red-shifted out- 
flow, which dominates the line profiles of CO. The lack of HCO + 
outflow wings as well as the absence of HCO + emission in most 
positions except around the central envelope, suggest that the 
swept-up gas is at a density of a few times 10 4 cm -3 or lower. 



3 RADEX is available online at 

http://www.sron.rug.nl/~vdtak/radex/radex.php 
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Outflow properties 





V " 

' max 


R 






M"- e 




Mem 




(km s _1 ) 


(AU) 


(M s ) 


(yr) 


(M G 


(Mo y^ 1 


(L s ) 












yr" 1 ) 


km s" 1 ) 




Red Lobe 


CO 2-1 


14.2 


5e4 


0.8/2.9" 


2.6e4 


9.1e-5 


9.8e-4 


1.2e-2 


CO 3-2 


14.0 


>2.5e4 


0.68 


>1.3e4 


4.2e-5 


4.5e-4 


5.4e-3 


CO 4-3' 


11.0 




>0.18 










CO 6-5 


10.5 


>2.5e4 


0.48 


2.2e4 


1.8e-5 


l.le-4 


8.0e-4 


CO 7-6 


12.1 


2.0e4 


0.4 


1.4e4 


2.2e-5 


1.7e-4 


1.6e-3 


Blue Lobe 


CO 2-1 


0.5 


1.2e4 


0.24 


1.2e4 


1.6e-5 


9.3e-5 


6.0e-4 


CO 3-2 


0.3 


1.6e4 


0.24 


1.5e4 


1.3e-5 


7.6e-5 


5.0e-4 


CO 4-3 


1.7 


l.le4 


0.19 


1.5e4 


l.le-5 


4.5e-5 


2.2e-4 


CO 6-5 


0.5 


l.le4 


7e-3 


l.le4 


5.4e-7 


3.0e-6 


2.0e-5 


CO 7-6 


1.1 


1.0e4 


le-2 


1.2e4 


8.6e-7 


4.3e-6 


2.4e-5 



" Actual LSR velocities; note that the quiescent gas is at Vlsr-5. 3 km s" 1 . SV max =abs(V ma x -V^R). Velocities are not corrected for inclination. 
* Corrected for inclination using the average correction factors of Cabrit & BertouJ j 19901) 
c Constant temperature of 100 K assumed for red outflow and 70 K for blue outflow. 
d Dynamical time scale : /j = R/SV mRX 

e Mass outflow rate : M = M/t^, not corrected for swept-up gas. 

f Outflow force : F co = MV^JR 

s Kinetic luminosity : L kin = 0.5M(5V max ) 3 /R 

0.8 M is derived for a temperature of 50 K. 2.9 M Q is derived for a temperature of 100 K 
' Outflow extends to far larger scales than map scale. 



(ratio=25). Fig.[9]shows that the line ratio increases for more ex- 
treme velocities, introducing a dependency of r w i ng on velocity. 
In addition, there is a large jump in r w j ng from 6 to 7 km s _1 , 
representing the transition between the optically thick quiescent 
and more optically thin shocked material. In the following anal- 
ysis, it is assumed that all shocked outflow emission (>7 km s , 
1.7 km s with respect to the systemic velocity of 5.3 km s _1 ) 
is optically thin for all transitions. The effects of the moderate 
optical depth of the lines are subsequently discussed. 

Fig. [10] shows the ratios of the CO 3-2/6-5 main beam an- 
tenna temperatures of line wings as functions of velocity for four 
different positions along the red outflow axis ((0,0), (-20,-20), (- 
30,-30) and (-40",-35")). The CO 6-5 data have not been binned 
to the larger CO 3-2 beam, so the comparison assumes similar 
volume filling factors of the shocked gas. Figure [TT1 shows the 
CO 6-5 spectra binned with a 10" beam and with a 20" beam. 
As can be seen, the differences between the spectra is negligible. 
Additional testing at other positions confirmed that the spectra 
do not differ by more than 20%. Ratios are only plotted if the 
emission in both wings is larger than 3cr. 

The kinetic temperature T^ n of the outflow can be derived 
by comparing the intensity ratios in the line wings from vari- 
ous tr ansitions with model line intensities of Ivan der Tak et all 
(2007). With the density assumed to be constant at a few times 
10 4 cirT 3 , the observed 3-2/6-5 ratios of 2-3 correspond to ki- 
netic temperatures of about 120 to 150 K. Outflow emission is 
slightly subthermally excited, especially the CO 6-5, with T ex 
ranging from 85 to 120 K. The rising ratios observed at the more 
extreme velocities in Fig.[10]correspond to lower kinetic temper- 
atures, but even the highest ratios of ~7 still indicate kinetic tem- 
peratures greater than 70 K. The variation of the optical depth 
with velocity, as seen in Fig. [9] could account for the rising ra- 
tios seen in Fig. [10] since a higher optical depth will result in 
a lower ratio for the same temperature. Even with the limit of 
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Fig. 9. The ratio of the main beam antenna temperatures of 12 CO 
with respect to 13 CO for the J =3-2 line at the (0,0) position. 
The ratios correspond to optical depths of 1.8 (ratio of 10) to 1.0 
(ratio of 25). 



The presence of a strong quiescent component at positions asso- 
ciated with the outflow is discussed in § 6. 

The wings of the isotopic lines provide an upper limit to the 
optical depth, T w ; ng , of the outflow. The ratio of Tmb of the 12 CO 
3-2/ 13 CO 3-2 line wings at the source position is shown in Fig. 
[9] as a function of velocity. For a constant density of 3xl0 4 cirT 3 
the observed ratios correspond to r w i ne of 1.8 (ratio=10) and 1 
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Fig. 10. The ratio of the 12 CO 3-2/ 12 CO 6-5 main beam antenna 
temperatures at various positions along the red outflow axis 
(solid:(Q,Q), dashed:(-2Q",-2Q"), dot-dashed: (-30", ,-30"), long 
dash:(-40" ,-35")) as functions of velocity. Average error bars at 
various velocities are given in the upper part of the figure at their 
respective velocity. They are applicable for all positions except 
(0,0), for which errors are a factor of 2 lower due to the longer 
integration time of CO 3-2 at the center position. 
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Fig. 11. The CO 6-5 spectrum at the 0,0 position rebinned at 10" 
beam (solid line) and with 20" (dotted line). 



T w ing=l-7, the infe rred kinetic tempe ratures will drop by only 
20% (see Fig. 4 of lJansen et al.l dl996l) ). 

Fig. PT2"lpresents the kinetic temperatures as functions of posi- 
tion along the red outflow for a constant density of 2x1 4 cm -3 . 
There is a clear trend towards lower temperatures at larger radii. 
Averaged between 7 and 9 km s _1 , temperatures drop from 170 
K close to the source to 80 K at a distance 40". However, the 
error bars, derived from the minimum and maximum ratio at ve- 
locities greater than 7 km s _I show that large variations are pos- 
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Fig. 12. The kinetic temperature of the outflow gas along the 
red outflow axis computed from the CO 3-2/6-5 line ratio for 
a constant density of 2xl0 4 cm' 3 . Diamonds are the tempera- 
tures derived from the velocity-averaged ratios. The limits are 
determined from the minimum and maximum ratios between 7 
and 1 1 km s _1 . The results are consistent with a constant temper- 
ature of ~ 100 K if variations in optical depth and density with 
distance are taken into account. 



sible. In addition to the optical depth effects discussed above, it 
is also possible that the higher ratios at more extreme velocities 
and larger distances correspond to lower densities n(H2). The 
low density of a few times 10 3 crrT 3 inferred for the surround- 
ing cloud raises the inferred kinetic temperatures from 80 K to 
100-120 K. Such a drop in density significantly would lower the 
T ex . The current observational data cannot distinguish between 
these two possibilities, but both options are consistent with a 
high kinetic outflow temperature of ~100 K that is constant 
with distance from the source. This temperature is si gnificantly 
higher than the 7\in and r ex of ~15 K assumed by lOlberg et al.l 
d 19921) for the CO 2-1/1-0 intensity ratios for the HH 46 outflow . 
Although some studies find T ex that low (Bachi iler et alj|200ll) . 
recent studies , at times including high-/ CO, also find hig her T ex 
and 7\i n (e.g. lHirano & Taniguchill2001llLee et alll2002l) . While 
a cooler outflow component is not excluded, our data clearly 
show the presence of warmer outflow gas. Kinetic temperatures 
as low as 50 K would require outflow densities in excess of 10 5 
cm -3 , which are excluded by the HCO + data. 

In a similar analysis, the blue-shifted outflow has kinetic 
temperatures of 70 to 100 K with a lower limit of 50 K, some- 
what cooler than the red lobe. In addition to a variation with ve- 
locity, the ratios also seem to vary with distance from the source. 
At larger distances, the CO 3-2/6-5 intensity ratio is almost a 
factor of 2 higher than at the source position. 

lHatchell et al. ( 1999) use a swept-up shell model to predict 
kinetic temperatures along the outflow axis and walls. Their 
predicted values of 50 to 100 K agree very well with our de- 
rived temperatures in both outflow lobes (Fig. [T2l . Such tem- 
peratures are much lower than cal culated for entrainment mod - 
els, which predict > 100 K (e.g. jLizano & Giovanardill 1995b . 
lArce & Goodman! (12002) compare different outflow models and 
show that at least some flows are best explained with the jet- 
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driven bow shock model. Fo r a thorough revi ew of outflow 
model s, see lArce et al.l (|2007). The models from lHatchell et all 
( 1999) predict an almost constant temperature along most of the 
outflow axis with increasing temperatures near the bow shock, 
the main site of energy deposition. This is consistent with our 
constant temperature along the outflow axis. 

5.2. Other outflow properties 

Additional properties of both outflow lobes are derived from 
the molecular em i ssion maps following the recipe outlined in 
Hogerheiide et al. (1998), in which the radii, masses, dynam- 
ical time scales, outflow force and kinetic luminosity are cal- 
culated (see Table |5j. The results are corrected for the inclina- 
tion of 35° found for this source dReipurth & Heathcotdfl99ll 
iMicono et al.l 1998]) using an averag e of the three correction fac- 
tors from ICabrit & Bertoutl dl990 ) (their Figures 5-7). These 
factors range from 1.2 to 2.5, especially for the mass, and are 
introduced to account for the difference in observed V max on the 
sky and the actual extreme velocities. Although excitation tem- 
perature variations are seen throughout the outflow that depend 
on velocity and distance from the source, we assume an average 
excitation temperature of 100 K for the red outflow lobe and 70 
K for the blue outflow lobe for the derivation of these parame- 
ters. 

The resulting values can be uncertain up to an order of mag- 
nitude due to the variations in the covered area and thus in radius, 
especially for the CO 4-3 (observed with FLASH). The mass es- 
timates of the low excitation CO lines (e.g., the CO 2-1) may be 
overestimated by up to a factor 4, due to the larger area covered 
and the assumption that the outflow is iso-thermal at 100 K. CO 
2-1 may be dominated by cooler gas with temperatures down 
to, say, 50 K, lowering the mass estimate to a lower limit 0.8 
M Q , only 25% higher than the masses found for the high-/ CO 
lines. Temperature differences are likely also responsible for the 
(smaller) difference in masses in the other lines. 

Even with the uncertainties in covered area, the dynamical 
time scales for these outflows all converge on 10,000 to 20,000 
years, with no difference between the red and blue outflow 
lobes. The values of Table HI are similar to the results found in 
ICabrit & Bertouil d 19921) and lHogerheiide et all d 19981) for other 
Class I outflows and also agree with the results from lOlberg et alJ 
(119921) who find L kin =4.5 xlO 3 L Q and dynamical time scales 
of4xl0 4 yr. Although our dynamical time scales are a factor of 
2 smaller, this can be accounted for by the smaller covered area 
in our observations. 

The main discrepancy with older studies is that our out- 
flow mass is up to an order of magnitude higher than that by 
lOlberg et al.l (1 19921). The origin of this difference is two-fold. 
First. lOlberg et am 19921) do not apply correction factors from 
ICabrit & Bertoutl (U990). Second, there is a large difference be- 
tween ass umed temperatures in the red lobe (15 K vs 100 K). 
Although lOlberg et alJ (11992) derive their temperature from the 
CO 1-0 and 2-1 emission, the presence and intensity of CO 4-3, 
6-5 and 7-6 outflow emission strongly constrain temperatures to 
our higher value of 100 K. As illustrated by the CO 2-1 example, 
this can introduce a significant difference in the masses. 

At fi rst sight, there is no diffe rence with the masses de- 
rived by IChernin & Massonl dl99ll). but they assumed a simi- 
lar temperature of 15-30 K as lOlberg et all (Il992l) . However, 
|Ch ernin & Massonl (1 19911) took an optical depth of ~5 for the 
outflow emission, which is not confirmed by our observed lack 
of outflow emission in the 13 CO lines. For our observed maxi- 



mum r w j ng = 1 .7, the outflow mass as derived from CO 3-2 would 
increase by a factor of 2 (x r/(l - e~ T )). 

From the dust map, a total (envelope + cloud) mass of 8 M 
is derived (see § 5.3.2). The total outflow mass in Table|5]can be 
as high as 3.2 M and thus consist of a significant portion (40%) 
of that total mass. If the low mass for CO 2-1 is adopted, this 
per centage drops to 10%. 

iBontemps et al.l d 19961) empirically derive a relation between 
the flow force, Fqo and bolometric luminosity 



log(Fco) = 0.91og(Lboi)-5.6 
and the envelope mass 
log(F co )= 1.1 log(M env )- 4.15 



(1) 



(2) 



Using the bolometric luminosity of 16 L and envelope mass 
of ~5 M , flow forces of between 3xl0~ 5 and 4xl0~ 4 M yr -1 
km s" 1 are expected. The parameters for the red outflow lobe 
indeed show flow forces for most CO transitions that agree with 
these predictions. Only the flow force of CO 2-1 is a factor of 2.5 
higher. If the lower mass of 0.8 M is used for this transition due 
to a lower outflow temperature of 50 K, the flow force similarly 
drops by a factor of 3.5, almost identical to the flow forces found 
for the CO 6-5 and 7-6. 

The observed flow forces for the blue outflow lobe are up to 
two orders of magnitude lower than those of the red lobe. This 
'missing' outflow force can be explained by assuming that the 
blue-shifted outflow escapes the envelope and surrounding cloud 
and thus interacts much less with the surrounding cloud due to a 
lack of cloud material in the path of this outflow lobe. 

6. Origin of the quiescent high-./ CO line emission 

6.1. Photon heating of cavity walls 

The analysis in § 4.2 shows that an envelope model derived from 
the dust emission is not able to fit the quiescent emission of the 
high-excitation lines by a factor of ~2.5, even though the low- 
excitation optically thin lines are well-fitted with standard abun- 
dances. The limits on the C I8 6-5, combined with the infor- 
mation that the I2 CO 6-5 emission from the envelope is opti- 
cally thick, implies that an additional heating component must 
be present that produces quiescent emission. There are several 
constraints on this component. First, it cannot originate in the 
inner regions of the protostellar envelope or be obscured by the 
envelope itself, because emission from such a component can- 
not escape through the optically thick outer envelope emission. 
Second, it must be (nearly) optically thin, since the observed 
C ls O 6-5 limit is already reached by the envelope itself. Thus, 
it cannot contribute significantly in mass. Third, it must be ex- 
tended since a similar warm quiescent component is also clearly 
seen at other positions covering the red-shifted outflow through 
narrow I2 CO 6-5 emission, e.g. at positions (-20" 10") and (- 
20", -20") (see Fig.|6]and[5). Observations of 13 CO 6-5 at other 
positions also confirm that this quiescent emission cannot be 
ver y optically thi c k (t < 0.4). 

ISpaans et al.l d 1995b investigated the influence of photon 
heating on the emission of high-/ CO lines, in order to explain 
the bright but narrow 12 CO and 13 CO 6-5 single spectra observa- 
tions of Class I sources (Hogerheii de et alU 998). In this process 
ultraviolet (UV) photons heat the gas of the outflow cavity walls 
to temperatures of a few hundred K, but are unable to dissoci- 
ate the CO molecules. The photons in this model originate in 
a ~ 10,000 K radiation field of the boundary layer in the accre- 
tion disk. It was proposed that dust present in the cavity scatters 
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the UV photons towards the cavity walls. The turbulent velocity 
of the gas at the cavity walls is low, thus explaining the narrow 
width of the CO 6-5 and 7-6 lines. This extra warm gas due to 
photon heating could account for the differences in the observed 
12 CO and 13 CO 6-5 emission on source and that modelled by 
just a protostellar envelope. 




be found within the outflow c avity of HH 46/47, th at likely orig- 
inate in the jet (see Fig. 3 in Velusamy et al. 2007). These posi- 
tions produce additional UV photons, creating a much more ex- 
tended PDR. The bow shock at the position of HH 47C, the two 
new jet/shock positions and the accretion shock boundary layer 
combined are excellent sources of the necessary UV photons. 
Figure Qj] shows a schematic overview of the red-shifted out- 
flow cavity and the proposed process of photon heating within 
the HH 46/47 outflow. § 6.2 discusses the path of the UV pho- 
tons in more detail, while § 6.3 presents the constraints on the 
UV field. 

6.2. Importance of the outflow geometry 

To illuminate a l arger area of the ca vity wall than possible by 
direct irradiation Spaan s~et al.l (1 19951) invoked scattering of the 
UV photons by the dust present in the outflow cavity, even 
though such scattering is dominated by anisotropic, mainly for- 
ward scattering. With a constant opening angle (see inset a of 
Fig- US, only a very small area of the envelope wall will be di- 
rectly impacted by UV photons. However, with densities of a few 
100 cirT 3 the mean free path of a UV photon until its first scat- 
tering event is ~65,000 AU, with the assump tion that on av erage 
an A v of 1 produces a single scattering event dDrainel2 003). This 
is significantly larger (~ 150" on the sky) than the area covered 
by the outflow cavity. 

A parabolic shape of the outflow cavity in the regions near 
the envelope, as seen in inset b of Fig. QjJ allows much more 
UV radiation from the accretion boundary layer to illuminate 
the cavity walls at larger distances of a few thousand AU. A 
more quantitative description of outflow shapes and its effect on 
the illumination of the outflow walls is beyond the scope of this 
paper. 



Fig. 13. Cartoon model of the HH 46 outflow on scales up to 
the bow shock (-60,000 AUx -60,000 AU), illustrating the pho- 
ton heating of the cavity walls that is responsible for the high-/ 
CO emission. UV photons are created both in the accretion disk 
boundary layer, as well as the bow shock and illuminate the cav- 
ity walls. Illustrative paths of the UV photons are shown with 
wigg ly lines. In i nset a , the method of photon heating proposed 
bv lSpaans et alj d 1995b using an outflow with constant opening 
angle is shown. There, dust must scatter UV photons to produce 
significant irradiation of the cavity walls. In inset b, a parabolic 
shape of the outflow cavity allows many more UV photons to 
directly impact the outflow cavity walls. Due to the expected 
long mean free path of UV photons within the outflow cavity 
(-65,000 AU), scenario b is assumed to be the more likely sce- 
nario. See the text for more information. 

Although Spaa ns et al.l (|T995) on ly include photons created 
in the accretion disk and thus limit their Photon Dominated 
Region (PDR) to a central region of up to 4,000 AU in size, this 
model can be extended if additional sources of UV photons are 
included. Such photons can be produced at bow shock positions 
of the outflow if /-shocks are present there dNeufeld & Dalgarnol 
119891) . Such photons can also heat the dust at the cavity walls, 
increasing the emission of the dust particles at continuum wave- 
lengths as presented in Fig. 2. The elongated continuum emis- 
sion of HH46 on a 1 ' scale in the south-east direction may indeed 
be caused by warm dust, instead o f more massive cold dust or 
CO emission. lVelusamv et al.l d2007l) show that apart from the 24 
fj.m hotspot of HH 47C, two additional 24 /im emission spots can 



6.3. Constraining the UV field 

The temperature of the cavity walls can be derived from the radi- 
ation field originating from these shocks. The radiation field Go 
at the cavity walls can be characterized by 

with G' the radiation field originating at each shock position i 
and f a geometric dilution factor to account for the difference 
in the UV emitting surface to the total illuminated surface of 
the cavity walls. A dissociative shock i produces 3xl0 7 (v s /100 
km s~ 1 ) 3 (no) photons cm" 2 s -1 as suming all energy is emitted 
at 10 eV (iNeufeld & Dalgarndll989l) . This is equivalent of a G 
of 30,000 assuming a density «o of 10 4 crrT 3 and a shock ve- 
locity v s of 220 km s _1 , derived for the bow shock dFernandesI 
2000). The average factor of /' over all shocks is assumed to be 
0.005, meaning that the working surface of each shock is -1/200 
in size compared to the cavity surface. However, not all shocks 
will have a velocity as high as 220 km s -1 . It is more likely 
that the velocities of the secondary shocks within the cavity are 
much lower. But similarly, the size of the working surfaces of 
the shocks, and thus the /' can be larger, resulting in a compara- 
ble overall strength Go, but of a different 'color' . With four spots, 
the cavity walls are illuminated by a total Go equaling -600 with 
an uncertainty of a factor of 2 to 3 due to the u ncertainty in the 
f. Using the results from Kaufm an et al.l d 19991 their Figure 1), 
the PDR surface temperature is then constrained to 250^1-00 K, 
sufficient to produce the observed lines. 
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[C I] 2-1 is detected near and at the bow shock position, 
but not at the outflow positions closer to the star, although the 
dynamic range in the [C I] data is small. This could indicate 
that the penetrating UV photons within the region of the out- 
flow lobe closer to the star are not able to dissociate CO sig- 
nifica ntly, constraining the color temperature of the radiation 
field. iNeufeld & Dalgarnol Jl989) show that shocks with veloci- 
ties less than 90 km s _1 do not produce CO dissociating photons. 
If all shocks have v s «80 km s _1 the estimated G' Q drops by a 
factor ~20. In that case, the total Go is not sufficient to heat the 
cavity walls to surface temperatures of 250-400 K, but only to 
about 100 K. Lower velocities are likely for the shocks observed 
insid e the cavity. Ho wever, the known shock velocity of 220 km 
s _1 dFernandesI l2000h is sufficient to produce CO dissociating 
photons, consistent with the observed narrow [C I] emission at 
this position. This situation is reminiscent of the observation of 
strong quiescent [C I] 1-0 emission in the supernova remnant IC 
443 ahead of the shock, ori ginating from p hotodissociation of 
CO in the pre-shocked gas ( Keene et al.lll996f) . The likely sce- 
nario for HH 46 is thus 

1. Non-dissociative UV photons are created in the boundary 
layer and secondary shocks, while the bow shock produces 
mainly CO dissociative UV photons 

2. Along the outflow axis, the cavity walls are heated to suffi- 
cient temperatures to produce the quiescent high-/ CO emis- 
sion. 

3. Closer to the bow shock, dissociation of CO becomes signif- 
icant in addition to the heating of the cavity walls, explaining 
the lack of CO 6-5 emission and presence of [C I] 2-1 emis- 
sion near and at the bow shock 

Observations of far-IR CO lines, such as are possible with 
the HIFI and PACS instruments on Herschel, are needed to con- 
strain the exact temperature of CO gas at the cavity walls. 

Slow (v s = 5-10 km s ) C-type shocks along the outflow 
cavit y walls are able to gener ate similar amounts of CO emis- 
sion dDraine & Rob erge 1984). However, the narrow nature of 
the line profile, as well as the presence of CO 6-5 emission over 
the entire area traced by the IR outflow, make this scenario less 
likely than the photon heating, although it could contribute some. 

7. Conclusions 

In this paper, we characterize the structure of protostellar 
envelope and the molecular outflow associated with HH 46 
IRS 1, as well as its immediate surrounding cloud material, 
through dust and molecular line maps. Broad and narrow CO 
lines are observed ranging in transitions from 2-1 to 7-6, 
including isotopologues. The three distinct components can be 
best described by the following model: 

- Envelope - The envelope of HH 46 with ~3-5 M Q (within 
T > 10 K) is one of the most massive ones found for a Class 
I low-mass protostar, but is densely concentrated toward the 
center (n oc r~ L8 ). The C ls O line emission from the envelope 
can be best fitted with a drop abundance of 2.7xl0~ 4 /i x 10~ 5 
above/below 30 K and below/above 10 5 crrT 3 . However, 
such abundanes are unable to reproduce the observed 12 CO 
and 13 CO 6-5 and 7-6 emission. The dense envelope 
itself is best traced by the HCO + 4-3 emission, which 
has very little outflow contribution and shows a spherical 
distribution. Densities in the inner few hundred AU of the 
envelope are high (>10 7 cm 3 ), with high optical depths 



of the HCO + 4-3 and all 12 CO lines. A C/H 2 abundance 
of a few times 10" 7 is found, which can be maintained by 
photodissociation of CO by cosmic ray induced UV photons. 

- Surrounding cloud - The surrounding cloud extends over 
more than 100" to the south-southwest but does not extend 
further than ~30" north of HH 46 IRS 1 where even no low 
excitation CO emission is found. Cloud conditions include a 
low density of a few times 10 3 cm -3 , derived from limits on 
the CO 6-5 emission at positions such as (30", -20"). The 
total column of CO is ~10 18 cirT 2 . 

- Outflow - The red-shifted molecular outflow, extending at 
least 40" from the source, produces strong molecular line 
wings up to CO 7=7-6 and heats the surrounding cloud and 
envelope significantly close to the star. Spatially, the red- 
shifted outflow lobe corresponds to the brig ht infrared out- 
flow lobe from Noriega-Cres po et"ai] {2004). Optical depth 
of the CO 3-2 outflow wing is less than 1.7. Kinetic tem- 
peratures of the red-shifted outflow are of order 100-150 K 
close to the star for flow densities of 2x 10 4 cirT 3 , but drop to 
80 K further from the central source if densities and optical 
depth remain constant. However, the data are also consis- 
tent with a constant kinetic temperature in the covered area 
if densities decrease to a few 10 3 at a distance >40" from the 
central source, as found for the surrounding cloud. 

Temperatures of both outflow lobes are significantly 
higher than the previously derived tempera ture of 15 K, 
but ag ree well with the model predictions of lHatchell et alj 
(1999) for a swept-up shell model. The high temperature 
causes the observed outflow mass to be significantly higher 
(almost an order of magni tude) than derived in older studies 
such as lOlberg et all d 1992b . Bright narrow [C I] is found 
near the bow shock, indicating that the bow shock produces 
CO dissociating photons. 

- Origin of high-/ CO - The emission seen in the higher ex- 
citation CO transition has three main origins. 

1 . The dense envelope produces optically thick emission in 
both CO 6-5 and 7-6, originating in the warm (T > 50 
K) inner envelope, accounting for roughly 1/3 of the ob- 
served line intensities on source. 

2. High-/ CO emission is detected in the red- and blue- 
shifted outflow wings at some positions along the out- 
flow axis. 

3. The bulk of the high-/ CO emission has narrow lines 
and is produced by photon heating. UV photons originat- 
ing in the bow shocks, jet shocks and accretion boundary 
layer heat the cavity walls up to a few hundred K. The 
lack of strong associated [C I] emission near the source 
indicates that the UV photons do not photodissociate CO, 
suggesting shock velocities lower than 90 km s such as 
could be present inside the cavity. CO dissociating pho- 
tons are limited to the region close to the bow shock. 

This paper shows that the addition of the high-/ CO emission 
lines as observed with CHAMP + provides new insights into the 
physical structure of the protostellar envelope and molecular out- 
flow of HH 46. The presence of narrow line emission in the 
high-/ CO lines throughout the outflow suggests that photon 
heating is an important process in HH 46. The high-/ isotopic 
CO lines, in particular those of C 18 6-5, have been key in con- 
straining the envelope model. The UV radiation implied by the 
high-J CO and [C I] observations should also have significant 
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consequences for the chemistry of other species and should en- 
hance radicals like CN and OH along the outflow walls. Future 
high-frequency observations using high=7 CHAMP + , Herschel 
and, in the long run, ALMA high frequency bands, will provide 
unique constraints on the interaction between outflow and enve- 
lope. 
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